Van der Waals interaction between an atom with spherical plasma shell 
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The van der Waals interaction energy of an atom with infinitely thin sphere with finite conductivity 
is investigated in the framework of the hydrodynamic approach. Thin sphere models the fullerene. 
We put the sphere into spherical cavity inside the infinite dielectric media, then calculate the energy 
of vacuum fluctuations in the context of the zeta-function approach. The interaction energy for a 
single atom is obtained from this expression in the limit of the rare media. The Casimir-Polder 
expression for an atom and plate is recovered in the limit of the infinite radius of the sphere. 
Assuming a finite radius of the sphere, the interaction energy of an atom falls down monotonic as 
I , third power of distance between atom and sphere for short distance and as seventh power for large 

I ■ distance from the sphere. Numerically, the interaction energy is obtained to be 3.8eV for hydrogen 

atom placed on the surface of the sphere with parameters of fullerene Ceo- We show also that the 
£\j i polarizability of fullerene is merely cube of its radius. 

■ PACS numbers: 73.22.-f, 34.50. Dy, 12.20.Ds 

S: 

Qs, • I. INTRODUCTION 

(N : 

The general theory of the van der Waals force was developed by Lifshits in Refs. [l|, 0| in the framework 
of statistical physics. In the case of interaction between particle and plate it is commonly referred to as the 
Casimir-Polder force [H . For small distance the potential of interaction is proportional to inverse third degree 
, of distance from the plate. For large distance the retardation of the interaction is taken into account and the 
potential falls down as fourth degree of distance. The last achievements in Casimir effect have been discussed 
in great depth in books and reviews 043- 

The van der Waals force is very important for interaction of graphene (graphite layers) with bodies @-[IB| 
and microparticles p^ - |20j . An understanding of the mechanisms of molecule-nanostructure interaction is of 
importance for the problem of hydrogen storage in carbon nanostructures [2l[ . The microscopic mechanisms 
' underlying the absorption phenomenon remain unclear (see, for example review [22|). 

In the present paper we use model of the fullerene in terms of the two dimensional free electron gas [23j 
which is usually called as hydrodynamical model. This model was applied and developed for the molecule 
Cqo in Refs. [24], [25|, for flat plasma sheet in Ref. (2(| and for spherical plasma surface in Ref. [27j • In 
the framework of this model the conductive surface is considered as infinitely thin shell with the specific wave 
number f2 = Anne 2 /mc 2 , where n is surface density of electrons and m is the electron mass. Since the surface 
is infinitely thin, the information about the properties of the surface is encoded in the boundary conditions on 
the conductive surface which are different for TE and TM modes. In the Ref. [27j it was shown that the energy 
of the vacuum electromagnetic fluctuations for surface shaped as sphere has a maximum for radius of sphere 
approximately equal to the specific wavelength of the model An = 2tt/CI. What this means is the Casimir force 
tries to enlarge sphere with radius larger then An and it tries to reduce the sphere with radius larger then An. 
The Boyer result [28| is recovered in the limit f2 — > oo. 

At the same time it is well known 29] that the energy of electrons in graphene has linear frequency dependence 
whereas in framework of the hydrodynamic model the energy of electrons is quadratic in the frequency. There 
is also another point that the electrons in the graphene have zero or very small effective mass. To describe 
correctly these unusual properties of electrons in graphene the Dirac fermion model was suggested in Ref. :30f . 
The electrons in this model arc described by (2 + \)D Dirac action with characteristic propagation velocity as 
Fermi velocity vp s» c/300 and very small mass gap m < O.leV. This model was applied for calculation of 
Casimir interaction energy between graphene plate and perfect conductor plane in Ref. [151 ] and recently in Ref. 
[20l ] for Casimir-Polder interaction energy between graphene and H, He* and Na atoms. 
It was shown that the Casimir energy for large distance between graphene plate and perfect conductor plane 
is decreasing by one power of the separation a faster than for ideal conductors, that is as (ma)~ 4 . If the 
mass of gap is zero at the beginning of calculations, m = 0, they obtained standard dependence a -3 . For the 
case of Casimir-Polder interaction energy between graphene and atoms [201 ] the hydrodynamic and the Dirac 
models give qualitatively different results. For the large separation the energy decreases with separation as 
a -4 which is a typical behavior of the atom-plate interaction at relativistic separations, but the coefficients are 
different. In the case of H, He* and Na atoms, the hydrodynamic model gives « 5 times larger coefficient than 
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the Dirac model. There is also interesting observation about mass gap parameter: the energy does not depend 
on the parameter for m < 10~ 3 / 2 eV and therefore the limit m — > is satisfied. 

There is another approach for van der Waals interaction based on the density-functional theory [3TI, HH and 
the local-density approximation [32| which has proved to be a very useful tool for calculating the ground-state 
properties of atoms, molecules, and solids. In framework of density-functional theory a number of studies of van 
der Waals interaction has been made fl3l [33l - l4ll j . The main problem in this theory is to find approximations 
for the exchange-correlation energy. The density-functional theory describes cohesion, bonds, structures, and 
other properties very well for dense molecules and materials. The theory fails to describe the interactions at 
sparse electron densities. The solution of this point by introducing the non-local correlations may be found in 
Refs. OH Ell- 
in the present paper the hydrodynamical model of fullerene is adopted - the infinitely thin sphere with radius 
R in vacuum and finite conductivity. To obtain the van der Waals interaction energy between an atom and this 
sphere we use the following approach which is due to Lifshits (see Refs. 0, H, [3 [l9| ) . We put the sphere inside 
the spherical vacuum cavity with radius L = R + d > R which is inside the dielectric media with coefficients 
fx, e. Then we find the zero-point energy of this system by using the zeta-function regularization approach, and 
take the limit of the rared media with e = 1 + AnNa + 0(N 2 ), where N — > is the volume density of the atoms 
and a is the polarizability of the unit atom. The interaction energy per unit atom which is situated d from the 
sphere is found by simple formula 

E n (s) = - hm 



n^o 4ttN(R + df 



where E(s) is the zeta-regularized energy with regularization parameter s. 

The paper is organized as follows. In Sec. [IT] we derive the boundary conditions for electromagnetic field on 
the infinitely thin conductive sphere as well as on the boundary of the cavity. Section IIIII is devoted to the 
construction of solutions satisfying the boundary conditions. The expression for the van der Waals energy is 
found in Sec. IIVI and it is analyzed in the limits of infinite radius of the sphere and for short and large distances 
between atom and sphere. Section [V] contains the numerical calculations of the interaction energy between 
hydrogen atom and the infinitely thin sphere with parameters of the fullerene Cqo . In the last section I VII we 
discuss results obtained. 



II. MAXWELL'S EQUATIONS AND MATCHING CONDITIONS 

Let us consider a conductive infinitely thin sphere with radius R in vacuum spherical cavity with radius 
L = R + d which is inside the dielectric media with parameters fi, e (see fig. [I]). We have two concentric spheres 




FIG. 1: The infinitely thin conductive sphere with radius R is located inside the vacuum spherical cavity with radius L = R + d 
into dielectric media with e, fi ^ 1. 

and we should consider the boundary conditions on two spherical boundaries. 

I. First of all let us consider a spherical boundary with radius L = R + d. Inside the sphere we have 
vacuum, e = fx = 1 and outside - the dielectric media with s, fj, 7^ 1. Assuming the spherical symmetry, the 
electromagnetic field is factorized for two independent polarizations usually called as te and tm modes. The 
Maxwell equations with oscillatory time dependence exp(— iuit) read 



= 0, divB = 0, 
= 0, divD = 0, 



(la) 
(lb) 
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where we should use the material equations D = e(w)E and B = /i(u;)H. For te mode to be obtained we 
express B from the first equation (fTa|) and substitute it into the second Eq. (Ilb[) 



B TE = --rotE TE , AE TE - — ^-/ieE TE = 0. (2) 
For tm mode to be obtained we express E from second equation (fTbf and substitute it into the first Eq. ([Taj) 



E™ = —rot B™ , AB™ - ^ M eB™ =0. (3) 

UJfJ,S c z 

We next expand solutions over spherical functions, Y/ m , and obtain the following expressions for te and tm 
polarizations, 

B^ = --rotE^, E^ = f(kr)LY lm , (4) 

UJ 

E™ = — rotB™, B™ = f(kr)LYi m , (5) 

where ck — Lo^fJIe. In the standard spherical vector basis, (e ri eg, e v ), we obtain in manifest form the modes we 
need 



= [Q,-^-,d 9 Yi m ,-ifd e Y lr , 



it 

sin ( 

r te ( c fui + -\\v c ( rf Y a v, ^iDLn v I 

a lm ~ H* + ±) ¥ lm, Oeilm, ; — T.Ou,Ilm , 

\uir tor Lursmt) 1 

B™ = [0,4^d v Y lm ,-ifdgY lm 
\ sin # 

e™ = -—(—Ki+mm,—m m ,-^- E d v Y lm ), 



where the function / obeys the following radial equation 



r+Ws'iff+i-w.o. ( 7 ) 



r 

The two independent solutions of this equation are the spherical Bessel functions ji{z) = 
^/ir/2zJ l+1/2 (z), yi{z) = ^Tr/2zY l+1/2 (z), where z = ruj^/eJI/c. 
At the boundary, r = L, the matching conditions read 

n • [B 2 Bi] L = 0, n • [D 2 - D X ] L = 0, (8a) 
n x [H 2 - Hji = 0, n x [E 2 - E^ = 0, (8b) 

where n = r/r is an unit normal to the sphere. We have to take into account also that k = oj/c inside the 
sphere r = L and k = ujy/JIe/c outside the sphere. The square brackets above denote the coincidence limit on 
the boundary r = L. 

II. The electromagnetic fields given infinitely thin conductive surface E in vacuum was considered by Fetter 
in Ref. (HI- The applications of this model for vacuum fluctuations of field see in Refs. [24T[27l | . The electrons 
of conductivity on the sphere produce current and the Maxwell equations read 

rotE-^H = 0, divH = 0, (9a) 

rot H + — E = 4ttJ, divE = 47rp, (9b) 
c 

where p = <5(x — xs)c, J = £(x — xs)j/c. Taking into account the equation of continuity and the Newton 
equations we obtain the following expressions for density and of charge and current on the boundary: 

? ? 

V r E», 3 = 1 E,|, (10) 



muj 2 mijj 



where the superscripts || indicate the vector components parallel to the surface E; e and m are the charge and 
mass of electron, and n is a surface density of charge. 
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As a consequence of the charge and current obtained above, the boundary conditions on the sphere with 
r = R read 

n • [H 2 - Hi] fl = 0, n • [E 2 - E^r = ^Vj • Ey, (11a) 
n x [H 2 - Hi] fl = - l —n x Ey, n x [E 2 - Ei] fl = 0, (lib) 

where k — lu/c and ft = Anne 2 /mc 2 is a specific wave number on the sphere. Because of the fact that the 
sphere is infinitely thin we may consider the Maxwell equations © in vacuum with zero right hand side and 
all information about sphere will be encoded in boundary conditions (jlip. An interesting treatment of this 
boundary condition is in Ref. [42| . 



III. THE SOLUTION OF THE MATCHING CONDITIONS 



Let us represent the radial function in the following way 

!fin = Oinji{kr), r<R 
font = a out ji{kr) + b out yi(kr), R < r < L (12) 
f e =a s hf\kr), r>L 

where ji,yi and are the spherical Bessel functions and k = uj/c inside the sphere, r < L and k = uj^/JjE/c 
outside the sphere for r > L. 

In this case the matching conditions (JU) and (TTT1) in manifest form read 

[rfout ~ rf in ] R = 0, 
[(rfout)' r - {rf in )' r - n{rfin)U = 0, 

[rfout - rf e ] L = 0, (13) 

MoutYr ~ -(rfeYrh = 0, 

for te mode, and 

[(rfout)' r ~ (rfin)' r ]R = 0, 

[(rfout) - (rfin) + -j^(rfin)r]R = 0, 

[rfout --rf £ ] L = 0, (14) 
M 

[(rfoutYr-— (rfe)'rU = 0, 
/is 

for tm mode. The solutions of these equations exist if and only if the following equations are satisfied 

' H(z e )9' TE --H , (z e )* rE = 0, (15a) 



fie \i 

' H(z e )y' TM + -H'(z e )y TM = 0, (15b) 



/fie e 

where z £ — z^JJie, z = kL = ujL/c; the prime is derivative with respect the argument, and 

*te(2) = J(z) + -J{x)[J{x)Y(z)-J(z)Y(x)\, (16a) 
k 

#tm(*) = J(z) + jJ'(x)[J'(x)Y(z)-J(z)Y'(x)}. (16b) 

Here J(x) = xji(x), Y(x) = xyi(x), H(x) — xh\ i(x) are the Riccati-Bessel functions, and x = kR. Therefore 
the functions we need (see next section) to obtain the spectrum of the energy read (we set \i = 1) 

Ete = H'(z £ )V TE - -^H(z £ )y' TE , (17a) 
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Etm = ff(z e )*' TM -4=#'(*0*TM. (17b) 



For e = 1, the result obtained in the Ref. [27| is recovered 



— TL = i\l-^J{x)H(x)\ =/7'tE I/-)- (1^0 



S TM = -i |l - ^J'(z)iT'(z)} = -if TM (k), (18b) 

for real value of k, and for imaginary axis k — > ik we obtain from above expressions the Jost functions in 
imaginary axis: 

S TE = i|l + ^s £ (x)ei(x)| =if TE {ik), (19a) 
Stm = -« jl - ^s!( a; ) e !( a; )| = -ifjM(ik), (19b) 
because H(ix) = (—i) l+1 ei(x), J(ix) = i l+1 si(x) and F(«x) = —i l si(x) — (—i) l ei(x), where 

/ 7TX / 2x 

sjfc) = J -jh+i/2{x), ei(x) = y —K l+1/2 {x) (20) 

are the Riccatti-Bessel spherical functions of the second kind. For the problem with z = to be avoided we 
multiply Ejm for z 2 

Ete = -i{H , (z e )9 T E--j=H(z e )*' TE y ( 21a ) 
Etm = -« 2 |i7(z e )^M-^^'(^)*TM|. (21b) 
On the imaginary axis k — > ik we obtain 

S TE = -^ e 'C s <0 $ TE-ej(* f! )*TE, (22a) 

Stm = z 2 |e ; (2 £ )$VM-^e;(z e )$TM|, (22b) 

$te = si(z) + —si{x)[si{z)ei{x)-si{x)ei{z)}, (22c) 
a; 

$tm = Sl {z) -Q-s'^x^izy^x) - s'^e^z)], (22d) 

where Q = fi-R, 2 = /ci, z e = z^/e, a; = kR and e = For e = 1 wc obtain 

Ste = heiik), Stm = z 2 / T M(«fc) (23) 

in accordance with Ref. [13] ■ 

IV. THE ENERGY 

Within the limits of approach suggested in Ref. [43| , the expressions for te and tm contributions in regularized 
zero-point energy read (w = kc) 

E TE (s) = _ ncC0S7rs ^sy v / dkk i-2s dklnJ ; TEi (24) 

E™(s) = - nccoS7TS / dfcfci-^lnSxM, (25) 
*■ 1=1 7o 

where the integrand functions are given by Eqs. (|22j) . The summations in these expressions begin with I = 1 
because for I = the electromagnetic modes (JSJ are zero. 
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The derivative of the regularized energy with respect to the distance d (E(s) = E TE (s) + E (s)) may be 
found by interchanging the derivative and summation with integration. In manifest form it reads 



d d E(s) = - 



He COS 7TS 



,2.s 



v />00 

Vi/ / dkk 1 - 



2s d, 



where 



Gte 
Qtm 



1 $' 



TE 



J TE 



Ve$TE e/(z £ ) ei(z £ )<S> TE ' 

$ tm _ 1 gKfs) 

*tTM 7ee;(z g ) 

■I'. 



TM e K*«) , 1/2 



*TM e '( z e) 2 2 V£ 



Let us consider now the rared media with e(iu>) = 1 + 4nNa{iuj) + 0(N 2 ), where a is polarizability of the 
atom and the density of the dielectric matter N — > 0. In this case the Casimir energy E(s) is expressed in terms 
the energy per unit atom E a (s) by relation 



/>oo 

E(s) = N E a (s)4n(R + rfdr + 0(N 2 ). 
Jd 

From this expression it follows that 

d d E{s) 



E a (s) 



lim 



JV^o 4nN(R + d) 2 ' 
and in manifest form we obtain the interaction energy per unit atom 

hCjJL 28 COS 7TS 



where 



G TE - 
Gjm = 



JTE 



w(R + d) 2 
fTE(ik) 



(26) 



(27) 



(28) 



e;(z)$TE e ( (z)$T E ' 



e ; ( z Wtm,z + e i ( z )*tm ^ e[ (z)<S>' TUz + e x (z)$ T M - 



With definitions of the functions $te and <I>tm we have the following relations 

$te = si(z)f TE (ik) - —sf(x)ei(z), 

fl 2 

$tm = si{z)f T u{ik) + —s\ (x)ei(z). 

Taking into consideration these expressions we express above formulas in slightly different form, 

Qs 2 (x)e 2 (z) 



G je = ei(z)si(z) 



Gtm = -e' l {z)s' l {z)-e l {z)s l {z) 1 ^ Q 



x hu{ik) 



s't(x)e' 2 (z) + S 't(x)e 2 (z) 



a, 



by separating the terms which have no dependence on the parameter Q = HR. By virtue of the fact that the 
Casimir energy is zero for an atom in vacuum (Q = 0) without boundaries, we subtract the terms with Q = 
and define the interaction energy by the following relation 



En = lim{E a (s)- lim E a (s)}. 

s— >0 Q— >0 

With this definition we integrate by part over k and arrive with the final formula (x = kR, z = kL) 



Eq. = 



hcfl 



dka(iuj) 



S 2 (x)e 2 (z) | s'i{x)e'i{z) + S 'j{x)e 2 {z)- 
fjE{ik) /tm(^) 



(29) 



(30) 
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where the Jost functions on the imaginary axes read 

f TE (ik) = l+-«,(a;)ei(a;), (31) 

f TM (ik) = 1 - ~*{(z)e{(z). (32) 

To perform computations one needs an expression for the atomic dynamic polarizabilities of hydrogen. In 
was shown in Ref. [44j that the precise expression for the atomic dynamic polarizability of hydrogen is given 
by the 10-oscillator formula 

«m = E ufk~a (33) 

k=l K ' a 

where Qh a are the oscillator strengths and uik,a are the eigenfrequencies. All these parameters may be found in 
Refs. (44||45j]. It was shown in Ref. [45| that the polarizabilities can be represented with sufficient precision in 
the framework of the single-oscillator model 

where a a (0) = 4.50 a.u. (1 a.u. = 1.482 • 10 _31 m 3 ) and ui a = 11.65eV for hydrogen atom. 

One can see from the expression p0[) that the energy is negative because the integrand is positive for arbitrary 
radius of the sphere, the wave number fi and arbitrary position of atom. The same observation was noted in 
Ref. [46] for ideal case. Let us consider different limits. 

1) In the limit of perfect conductivity, il — ¥ oo, which we call the Boyer limit, we obtain 

»(R + ci) 2 ^ J a '|s,(x)e,(z) si^MM J 

2) The limit of infinite radius of sphere, R —} oo, with fixed distance, d, between the surface of sphere and 
an atom requires more machinery. One cannot merely interchange the limit and summation and integration in 
above expressions (I3U1) and ([33]) because in this case the integrand has no dependence on the I and the series is 
divergent. Indeed, in the limit of infinite radius of sphere 

2si(x)ei(z)\ R ^ co = +e~ kd , 2si(x)e l (x)\ R ^ 00 = +1, 
2s' l {x)e' l {z)\ R ^ 00 = -e- kd , 2 S ' l (x)e' l (x)\ R ^ 00 = -1, 
2s' l (x)e l (z)\ R ^ 00 = +er kd , 

and the sum over I is divergent, 

-2kd 

2 n (R + < Q*£ 



E ^ = ~ oZ71T~T' ^ 2 ^ y / dka(iw) — R -> oo. (36) 

2fc 



In order to obtain the correct expression for the energy in the limit R — > oo we change the variable of 
integration k — > vk in Eqs. (|3"0")) and (1331) 

hen " 2 J gfMefM 4(^)et(^) + ^Me?(^)^ v 

^ n = cp_l^2 2^ " / dka ( iuu ) \ f el \ + 7 — rr—\ — f ' ( 37 

ir(R + d) \ l Jo i fTE(tkis) jTM^kiy) 

E B = hc fr 3 rkdkafkju) [ S 1M^M ^y 2 A^) + ^x)ej{ V z)^ ^ 

7r(i? + d) 2 ^ 7 I si(i>x)ei{vx) s' l (vx)e' l (vx) 

and use the uniform expansion for Bessel functions (see Ref. [13|). We obtain the following expressions 



En = - W r dba^e-MiW-iW] (^M + 1 + « (*) ■ „ 1 (39) 

7r(i? + d) 2 ^ J y ' I 4m Apxzt{x)t(z) J 

, OO 



^ 3 r dAAa^Je-^W-"^ ( ^ + + •••}> (40) 



?r(i? + d) 2 ^ 7 v ; [ 2 2zt(z) 
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where p = 1 + 2v £ t (x) > w = 1 + 2 ^f l > = l N 1 + x2 > n( x ) = Vl + x 2 + In 1+ J 1+x5 and x = kR, z = kL = 
k(R + d). In the limit of R — > oo, the integrands in above both expressions have the same form and the main 
contribution to the energy comes from the first term of uniform expansion, 

F= v hc 9 2 y 3 r dyy e^M^iv)] 

rSL nc 2 (R + d) 2 ^ J y 2 u 2 +q 2 ut(u) ' 1 ' 

where u = y(l + d/R), q a = k a R and we changed variable k — » y = kR. Here the single-oscillator model for 
polarizability (|34p was taken into account. 

Next, the sum over I we represent in the following integral 



3„-2vS 



1 f°° 27 + Ue-^+V + 5e~ 4 ( t+5 ) - e -6(*+ 5 ) . 2q a t 




,_ 1 y 2 v 2 + q 2 a ~ ^hy Jo e 3( t +*)( e - 2 ( t + 4 ) - 1) 4 sm y [ ' 

Assuming this expression we interchange the limit R — > oo and integrals over y and t and obtain 
where 

1 r°° ( 1 i + + ^ 

(44) 

and v — dk a . Let us consider large distance, d, between the plate (sphere of infinite radius) and an atom, 
dk a S> 1. In the limit of v —> oo we obtain that S = 1 and therefore the Casimir-Polder (~ d~ 4 ) energy, 

is recovered. For small distances, dk a <C 1, we change the variable t — > t = t/2v and take the limit of v — > 0. 
In this case we obtain that S = irv/3 and the energy has the form ~ d~ 3 , 

E = -^§^, (46) 
as should be the case. The plot of the S as function of variable v = dk a is shown in Fig. [5] 
S 



0.6 



0.4 



0.2 

dk fl 

2 4 6 8 10 

FIG. 2: The plot of S as the function of the v = kd a . It tends to unity for large v (E ~ d~ 4 ) and it is linear over v (E ~ for 
small distances between an atom and plate. The relation of the energy and S is given by Eq. II43H . 

3) Let us analyze the energy for large (rf> k~ , d 3> i?) and small (c? <C k^ 1 , d <C R) distances between the 
sphere and an atom for finite f2 and R. In the case of large distance, d — > oo, of an atom from the shell we use 
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Eq. (|30|) . We change integrand variable k = y/d, next take limit d—too, and then we take the integral over y. 
The main contribution comes from the first term with 1 = 1: 

Su-Zl^ + Sml (47b) 



d 3 [3(3 + Q) 3 

^ (a) = W — vr^ — e dy ' (47c) 

where a 2 = Qd 2 / R 2 = d 2 Q/R. The first term in above expression (|47h>|) comes from TE mode and second - 
from TM polarization. The function F increases monotonically from zero for small a (d 2 <C R/ft) to unity for 
large a (d 2 ^> R/ft). In the case of a <C 1 the function F(a) ~ 27r v / 6a/23. Therefore, in the limit of ft — > 0, the 

energy Eq — > as should be the case. 

Assuming a finite conductivity, SI 7^ 0, and large distance d 3> fc^ 1 , <i R,d^$> yj R/Q we obtain that 

i? 3 f 46 

<3n 



d 3 (3(3 + Q) 3 
and we arrive with expression 

(48) 

Taking into account the Casimir- Polder interaction energy of two atoms with polarizations a and a / , 

23 ^(0)^(0) (49) 

we observe that the sphere with finite conductivity has static polarizability 

53Q+138. 

at = — - R 6 . (50) 

1 46Q+138 v ; 

To analyze the energy for small distances we use the following representation for the series 

,2 „-1vb 



^^ 2 + g 2 l + £ 4(Q 2 + a 2 y2)j / o |/ + 2g / ^y+g 7 &m y J 



(51) 



where /(x) = e 3 (*+ x )/(l — e 2 ( s + x )y The first and second terms give the d 3 contribution and the last term 
gives contribution ~ d. Taking into account these expressions we obtain that 

hca(0)k a 

E = MT~ ( 52 ) 

as should be the case, because close to the sphere we observe flat surface. 

V. NUMERIC ALS 

We analyze numerically the following expression for the energy (x — kR, z = k(R + d)) 

' s 2 (x)e 2 (z) , s' 2 (x)e' 2 (z) + s' 2 (x)e^(z) v ^\ 



z^-^Tdy^i dka ^ 

where the Jost functions in imaginary axes read 



/te(«&) frM(ik) 



(53) 



f TE (ik) = l + ^ Sl (x) ei (x), (54) 
f TM (ik) = 1 - js^xy^x), (55) 
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and polarizability of atom has the single-oscillatory form 

«M = -4^- (56) 

In the Boyer limit D, — > oo we obtain 



E F 



ir{R + df ^ 



V:L dkkai ^\^m wm r (57) 



For simplicity we extract as a factor the Casimir-Polder expression for the interaction energy of an atom with 
plate, 

3hca(0) 

&(l,B — 87rd 4 n ' B ' ' 

and we will numerically calculate the dimensionless quantities 

3(l + r)2^^7 y 2 + ^\ / TE ( % ) + /tm (iy) J' 1 j 

B m + r) 3 ^ Jo y 2 + q 2 a\sL(y)ei(y) <(yK(y) J' 1 ' 

where z = (1 + r)j/, <7 a = u a R/c, r = d/R, y — kR and 

/ TE fe) = l + -si(y)ei(y), (61) 
2/ 

/tm (*y) = ( 62 ) 

We use l/k a as the unit of measurement of length and therefore the function S depends on the three param- 
eters: £l/k a , q a — Rk a and dk a . The numerical analysis of the function S for Q/k a — 2.44 • 10~ 2 (molecule Cqo) 
and VL/k a = 1 is shown in Fig. [3] 




FIG. 3: The plot of S as the function of the dk a e (0, 2) and Rk a € (0.02, 2) for a) Cl/k a = 2.44 • 10" 2 and 6) Ct/k a = 1. 



Let us consider the interaction energy between hydrogen atom and molecule Cgo- For this molecule we have 
[H: R = 3.42A = 0.342nm, Q = ClR = 4.94 • 10" 4 and fl/k a = 2.44 ■ 10" 2 . The polarizability of hydrogen 
atom within the single-oscillator model reads 0, 0, El] a a (0) = 4.50 a.u. (1 a.u. = 1.482 • 10~ 31 to 3 ) and 
0J a = Il.tteV = 17.698 • l0 15 Hz (k a = 0.059nm _1 , A a = 106.4nm) where cj/c = k = 2tt/A. Therefore, 
q a = k a R = 0.0202. 

Taking into consideration all the numerical values of parameters we represent the energy for this system in 
the following form 



En{eV) = -Sn{q a ,r), 

d 4 (nm) 



(63) 



11 




FIG. 4: The plot of S as the function of the distance d between an atom and the sphere. Thin curve is the energy for the case 
R — > oo (Casimir- Polder energy for plate), middle thickness curve is the case of the molecule Ceo, an d the thick curve is the case 
of ideal sphere (f2 — > co). In the figure b we compare the energy for the plane with the energy in the sphere case. 



where the energy is measured in eV and the distance is measured in nanometres. The numerical simulations 
for the function S are shown in Fig. @] and the energy Eq in Fig. [5J The radius of the hydrogen atom is 
rjf = 0.053nm. For this minimal distance, d = r#, we have numerically E = 3.8eV. In the case of plate with 
hydrogen atom we obtain 6.4eV . In the interval of distances from the hydrogen atom radius rn up to 5r# the 
energy is approximated by the following expression 



E n (eV) 



The same dependence was observed in Ref. [2C 
For large distances we obtain from Eq. (|48|) 



En{eV) 



0.00013 

d 7 / 2 (nm) 



0.0095 



(64) 



(65) 



This expression approximates the exact one with error 10% starting with distance d — 50nm. The Eq. (|50[) gives 
the static palarizability of the fullerene a p (0) — R 3 = 4 ■ 10~ 29 m 3 . This expression is close to that calculated 
in Ref. [49( where the authors obtained a p (0) — 7 ■ 10" 



" 29 m 3 . 



£n(eV) 




-rrrr 



d (nm) 



a) 



£ I2 (eV) 



-0.006 




1.0 



— u d (nm) 
2.0 



b) 



FIG. 5: The plot of the energy Eq as the function of the distance d between the sphere and the hydrogen atom. In the figure a) 
we show the energy starting from the distance d = 0.053(nm) (the radius of the hydrogen atom). In the figure b) the energy in 
large interval is shown. 



VI. CONCLUSION 



In the foregoing, we have obtained the analytic expression for the Casimir-Polder (van der Waals) energy for a 
system which contains an atom or microparticle and infinitely thin sphere with finite conductivity which models 
a fullerene. We used the zeta-regularization approach and for renormalization we used a simple physically 
reasonable condition - the energy should be zero for an atom alone without a sphere. The conductive sphere 
with radius R is characterized by the only parameter ft — Aitne 2 / mc 2 with dimension of wave number, where 
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n is the surface density of electrons. The limit Q — >■ oo corresponds to the ideal case considered by Boyer |28j |. 
The microparticle is characterized by the only parameter, polarizability a. 

The expression obtained reproduces in the limit R — > oo the Casimir-Polder result for an atom and plate 
(see Eqs. V For small distances we have d" 3 dependence and far from the plate we obtain d~ 4 due 

to retardation. For finite radius of the sphere we have different behavior of the energy. Close to the sphere, 
d -C l/k a and d <C R, we have the same d~ 3 dependence as in the Casimir-Polder case and far from the sphere 
we obtained d~ 7 dependence given in Eq. (l4"8"l) . This expression is valid for d ^> l/k a and d ^> R. For the 
interval Th < d < 5r#, where m is the radius of the hydrogen atom, the energy is approximated by d~' 7 1 2 
dependence. We also note that the finite conductivity decreases the energy in comparison with Boyer case which 
may be observed in Fig. 2J 

Application to the molecule Cgo with hydrogen atom is plotted in Fig. [5] For closest distance atom from the 
fullerene, which is radius of hydrogen atom rjj, the energy is 2>.&eV which is two times smaller then for the case 
of hydrogen atom with plate. Away from the fullerene (in fact larger then 50nm) the energy falls down as d~ 7 
(see Eq. (|65|) ) which is in three orders of magnitude faster then for the Casimir-Polder case. This dependence 
corresponds to the Casimir-Polder interaction atoms for large distance. Taking into account this analogy we 
obtain the polarizability of fullerene (Q = f2i? = 4.94 ■ 1CT 4 < 1) 

= 53Q+ 138 3 ^ 3 = _ 10 _ 2 9 m 3^ 
; 46Q + 138 

This expression is close to that calculated in Ref. [4!| where the authors obtained a p (0) = 7 • 10~ 29 m 3 . 

In the paper we considered the interaction energy in the framework of the hydrodynamical model. As it was 
noted in the Introduction this model does not describe correctly graphene and therefore the systems made of 
them such as fullerenes. The model which describes graphene more precisely is the Dirac model. Nevertheless, 
using the calculations within the Dirac model which was made in Ref. [2(| as the base we expect that the 
interaction energy in framework of the Dirac model will be in five times smaller at large distance between 
fullerene and an atom. The dependence on the energy for large and small distances between fullerene and an 
atom is expected to be the same. 

There is another question which was not considered in the paper but which is very important for condensed 
matter physics. It is interesting to obtain the adsorption energy of the hydrogen on the Cqo at the physical 
equilibrium distance. This question is very important for the problem of storage of hydrogen in carbon nano- 
systems (see Ref. [2l| and review [22J). We plan to investigate these questions in the future works. 
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